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I.  INTRODUCTION 


A  field  of  primary  interest  in  Operations  Research, 
especially  as  applied  to  supply  support  problems,  is  the 
study  of  inventory  models.  Inventory  theory  is  the  appli¬ 
cation  of  mathematical  optimization  techniques  to  the  problem 
of  deciding  when  and  in  what  quantities  certain  goods  should 
be  purchased  to  meet  the  demand  of  customers  while  minimizing 
time-average  cost  [1]. 

The  fundamental  inventory  model,  usually  called  the 
economic  order  quantity  (EOQ)  or  Wilson  model,  presumes  that 
all  costs  and  demand  are  fixed  over  time,  that  resupply 
occurs  after  some  deterministic  leadtime  following  placement 
of  an  order,  and  that  demand  will  be  met  without  incurring 
shortages  of  stock,  all  of  this  for  a  single  activity  hold¬ 
ing  only  one  good  or  commodity  (see  the  sketch  below) . 


_\l 

Purchase  Q. 

T"  u^its  o-P 

f 

r*wen4-*v*  M 
He\d 

\ 
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More  complex  and  realistic  models  have  been  developed  by 
relaxing  some  or  all  of  these  assumptions  [1] .  Nevertheless, 
the  majority  of  the  problems  that  have  been  solved  concern 
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finding  optimal  policies  for  single  activities  since  multi¬ 
echelon  problems  are  more  difficult  to  formulate  and  solve. 

For  many  practical  applications  inventory  models  for  a 
single  activity  suffice.  But  for  some  large  retail  firms 
and  military  supply  systems  which  have  centralized  procure¬ 
ment  for  a  system  of  retail  outlets  or  stock-points,  a 
mathematical  model  which  accounts  for  the  mult i- activity , 
multi-echelon  nature  of  the  inventory  system  is  more  appro¬ 
priate.  A  classic  principle  of  optimization  theory  is  that 
any  allocation  of  resources  for  a  multi-component  system 
based  on  optimal  solutions  for  each  component  operating 
independently  can  be  no  better  than  the  solution  found  for 
the  system  as  a  whole  by  subordinating  the  roles  of  the 
independent  components.  It  is  the  author's  belief  that  most 
mult i- echelon  inventory  systems  tend  to  operate  sub- op timally , 
because  the  above-mentioned  independent  component  approach 
is  used  to  solve  for  stocking  policies. 

This  paper  presents  a  simple,  deterministic,  multi¬ 
echelon  model  which  is  optimized  on  a  "system"  basis.  The 
model  is  an  extension  of  the  EOQ  with  two  echelons,  as 
depicted  below. 
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The  upper  echelon  consists  of  a  single  activity  which  could 
represent  a  centralized  procurement  and  warehousing  facility. 
The  lower  echelon  has  m  independent  retail  outlets  or  stock 
points  which  replenish  from  the  upper  echelon  and  meet  the 
demands  of  the  customers  outside  the  system.  A  minimum 
time-average  cost  solution  for  the  system,  presuming  no 
stockouts  are  allowed  and  subject  to  the  system  cycle  length 
being  an  integer  multiple  of  the  cycle  length  of  each  lower 
echelon  activity  (possibly  a  different  multiplier  for  each 
activity),  is  derived  in  Chapter  II.  In  Chapter  III,  an 
algorithm  is  developed,  to  find  the  best  choice  of  integer 
multipliers.  The  case  in  which  customers'  demands  may  be 
backordered  is  taken  up  in  Chapter  V.  Some  examples  of  the 
no-stockouts-allowed  model  solution  are  presented  in  Chapter 
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IV  with  an  analysis  of  the  cost  benefit  available  in  the 
use  of  this  model  over  the  use  of  independent  EOQ  models 
for  the  activities. 

The  work  in  this  paper  is  basically  an  extension  of 
the  two-echelon  model  with  one  activity  at  each  echelon 
sketched  in  [2]  and  presented  in  greater  depth  with  the 
finite  production  rate  and  backorders  cases  in  [3] . 

Several  authors  have  attacked  special  cases  of  the  multi¬ 
echelon  problem.  A  probabilistic  model  for  repairables  is 
presented  in  [4] .  A  probabilistic  two-echelon  model  is 
solved  in  [5]  using  a  dynamic  programming  technique. 

The  deterministic  models  in  [2]  and  [3] ,  while  they  are 
fine  in  themselves,  solve  problems  of  completely  specified 
size.  The  stochastic  models  of  [3]  and  [4]  are  not  multi¬ 
echelon  extensions  of  the  probability  models  presented  in 
Chapter  Four  of  [1] .  This  paper  is  a  natural  two-echelon 
extension  of  inventory  theory's  fundamental  model  with  an 
arbitrary  number  of  activities  at  the  lower  echelon.  It 
is  hoped  that  it  also  provides  a  starting  point  for  further 
development  of  more  complex  deterministic  models  and  more 
general  stochastic  models  such  as  extensions  of  those  in 
Chapter  Four  of  [1] . 

The  author  presumes  in  his  presentation  that  the  reader 
is  familiar  with  the  EOQ  model  and  its  backorders- allowed 
extension. 
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II.  THE  NO-STOCKOUTS -ALLOWED  MODEL 


A.  MODEL  FORMULATION 

Suppose  that  there  is  an  inventory  system  consisting  of 
two  echelons  of  inventory,  performing  the  procurement,  the 
holding,  and  distribution  of  a  single  good.  The  upper  echelon 
is  a  single  activity  which  orders  and  buys  the  good  from  an 
outside  source  and  holds  it  in  inventory  for  further  distri¬ 
bution  within  the  system.  The  lower  echelon  has  m  independent 
activities  which  order  and  resupply  from  the  upper  echelon 
and  hold  inventory  to  meet  the  demands  of  their  customers. 

In  setting  up  the  cost  equation  for  the  system  the  following 
are  assumed: 

1)  All  orders  are  filled  immediately  (a  deterministic 
lead  time  has  no  effect  on  the  optimal  policy) ; 

2)  all  demands  are  met  immediately; 

3)  the  demand  rate  is  constant  and  continuous  over  time 
but  may  be  different  for  each  lower  echelon  activity; 

4)  the  goods  are  purchased  only  by  the  upper  echelon; 
and  the  price  is  constant  over  time; 

5)  ordering  costs  are  constant  over  time; 

6)  inventory  holding  costs  are  products  of  the  annual 
interest  rate  and  the  purchase  price  of  the  good 
(following  [1]);  and 

7)  the  cycle  length  (time  between  orders)  will  be  constant 
for  each  activity;  and  system  cycle  length  will  be  a 
positive  integer  multiple  of  each  lower  echelon  activ¬ 
ity's  cycle  length;  the  number  of  complete  sub-cycles 
for  each  lower  echelon  activity  may  be  different. 
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The  sketch  below  is  a  schematic  of  the  material  flow  in 


this  inventory  system.  The  construction  of  the  mathematical 
model  of  the  system  parallels  that  of  the  EOQ  model. 


Throughout  this  paper  whenever  a  summation  sign  ( S  )  is 
used,  it  is  assumed  that  sum  is  taken  for  the  index  variable 
ranging  from  one  to  m  unless  it  is  specifically  stated 
otherwise . 

Activity  0  is  the  upper  echelon  activity.  Activities  one 
through  m  are  the  lower  echelon  activities.  The  subscript  of 
a  variable  identifies  the  parameter  or  variable  with  its 
associated  activity.  The  A^,  1^,  and  are,  respectively, 
the  ordering  cost,  interest  rate,  and  consumer  demand  rate 

Vi 

for  the  1  activity.  These,  along  with  C,  the  unit  price 
of  the  good,  are  the  system  parameters.  The  Q^,  T^,  and  n^ 
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are,  respectively,  the  ordering  quantity,  the  cycle  length, 
and  the  number  of  complete  cycles  over  Tq  for  the  il  activity; 

and  are  continuous  variables  while  n^  is  a  positive 
integer  variable. 

Let  Z.  be  the  total  costs  of  running  the  inventory  system 

for  one  system-cycle,  T  .  Since  each  cycle  is  comprised  of 

th 

n^  sub-cycles  for  the  i  activity,  total  ordering  costs  are 
Aq  +  Z.  "i At  .  Total  purchasing  cost  is  merely  CQq,  for  the 
goods  are  bought  only  at  the  upper  echelon.  Let  B^(t)  be 
the  quantity  of  inventory  on  hand  at  time  t.  Then  the  total 

VY>  I 

holding  costs  are  IiC  ]  Bj.ttlot.  .  Putting  the  parts 

together ; 

T0 

V  =  A0  *  ♦  CQo  +  Is  TiC  \  BiCtJdt. 

^  it  O  c 

A  typical  realization  of  the  on-hand  inventories  within  the 
system  over  one  cycle  is  shown  in  the  following  sketch. 
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The  assumptions  of  the  model  result  in: 


T  =  n.  T.  for  all  i; 
oil 

for  all  i;  and 


Q0  *  £  ni  Qi  • 

.i 

Now  consider  IiCJBiW  Jt 


(2-1) 


for  l*c£vv^.  Since 


the  n^  sub-cycles  are  identical  and  the  area  under  the  on- 

hand  curve  for  one  sub-cycle  is  just  that  under  a  triangle, 

T: 

a  Tv 


Lc/b-CtUt  ^  nJ-.C  y BiWdt*  nvrtC 


For  the  upper  echelon  activity  B0  (t)  can  be  decomposed 
into  the  sum  of  the  inventories  held  for  each  of  the  m  lower 
echelon  activities.  If  B*  (t)  is  the  inventory  held  for  the 
it^1  echelon,  then  Bq  (t)  =  E  B^  (t)  where 


B>)=  {n,-i  -  Intf-rJ 


Q; 


and  ts 


the  greatest  integer  less  than  or  equal  to  X. 
Then 


--  J  (n-.-l-rnt[%])Q-.dt 


i'O 

Summing  over  i  gives 

To 


J Lt'id  t  =  L  z  w-Cnt-U. 


2  = 


It  follows  that  the  total  costs  per  cycle  are 
A0+  +  CC±°  +  H,  i  IiCGii  +  L  •£ 
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B.  MINIMIZATION  WITH  RESPECT  TO  Q 

To  compute  the  time-average  cost  for  the  system,  divide 
the  cost  per  cycle  by  the  cycle  length,  TQ  =  T^.  Thus, 

K--  +  I1I1C& v  Z  (.1-2) 


Although  this  is  a  function  in  2m  +  1  continuous  vari¬ 
ables,  it  can  be  reduced  to  a  function  in  one  continuous 
variable  by  using  the  equations  (2-1)  that  arose  from  the 
model's  assumptions.  Since  Q-  =  D •  T  •  and  T  =  n-  T.  = 

r  111  Oil 

n^  Tj  for  i,j  =  l,...,m;  equation  (2-2)  becomes 


K- 


A»D,  T  AiO, 

n.<5,  +  L  *.<3, 


C  L£i.  + 


t  IlCDl^.Q,  y  T>Cn,Q,(n--l)  D; 
L  ar>,ru  L  2  D,m 


(2-3) 


In  this  form,  K  will  be  referred  to  as  K(Q,  ,  n,  , .  .  .  ,  nm)  . 

For  a  given  choice  of  n, , . . . ,n  K(Q,)  will  represent  equation 
(2-3). 


Taking  the  derivative  of  K(Q,)  gives 

fill  r  .  AeD,  _  y  Ay.  fly  D,  y  T;C  PLn,  7  I.C,  n.oyQD,. 

dQ.  rv.Q,1  **  n.Q,^  +  ^  2D,|^II  ¥  ^  2D,  fly  ,  (2-4) 


and 


Now , 


JK 

d<3. 


O  implies 


"  2  D,4  { A*  +[nL  At  j 

1 

 l/^l.Uy-DDy 

*hi>i  1 

'  l  J 

A^k  s 


240D, 
n,  Q,3 


X^Aj.  > 


O 


(z-S) 


since  all  of  the  parameters  and  variables  are  positive. 

A 

Therefore,  Q,  gives  the  minimum  of  K(Q,). 

A 

Substituting  Q,  given  by  equation  (2-5)  into  equation 
(2-3)  gives 
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K(Q-)= 


Z 


I.K-DPi  +  TvPl 


'  V* 

'  +  C  I  Pi. 


(I'O 


Let  (Q,  ,  A,  ,...  ,  nm  )  be  an  optimal  value 
of  (Q.jO., . ,  a„,  )  .  Then 

Theorem  1-1:  K(Q« >n>)  —  > )  =  K (Q» »  n'  > **  *  J r'M 
Proof:  Since  K(Q*  *'*•••> -  rv‘‘v'  ^  (&»»*•  } 

K  (Q^a*,...,  KCQ.,  Now,  for  each  fixed 


KC Qon,,...,nm)<  K(Q., *.,—>*»«»)  for  all  Ql. 

Then  K  (Q*>¥'*>" ‘  >  ^  )  £  K  (Q»>nr>— >Am  )  for  all  Q^. 

Thus,  K  (.Q*A*>  — >n"  )  1  KCOr,.~,0.  Finally, 
kcqV^O  *  K(Q.,or,...,o^) <  k(q*a*...,iO  which 

implies  the  result. 

y\ 

Corollary  1-2:  Q  ,*  =  Q»  •  •  •  >  )  . 

C.  CLOSING  REMARKS 

*  *  * 

Given  Qt  ,n  ,...,nm  the  remaining  components  of  the 
optimal  operating  doctrine  can  be  derived.  Let  starred 
variables  represent  the  optimal  values  of  the  variables. 
Then 


T*  « 

5,  <3* 

• 

) 

(2-7) 

ry  = 

A,  -r~  * 

Al  1  * 

•pftr  <X.  1  1  L  j 

(2-8) 

* 

Qi 

Dili* 

for  ali  l  j 

(2-9) 

r/= 

*j  A  ni 

(2-10) 
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a *  --  i  <  qc  . 


(2-11) 


In  this  chapter,  the  no-stockouts -allowed  model  has  been 
developed,  and  the  values  of  the  continuous  variables  which 
minimize  the  time- average  cost  function  (given  any  set 

of  the  integer  variables)  have  been  found. 

In  the  next  chapter,  a  method  for  finding  the  optimal  values 
of  the  integer  variables  is  derived  which  will,  by  using 
equations  (2-7)  through  (2-11),  give  the  optimal  policy  for 
operating  the  inventory  system. 
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III.  FINDING  OPTIMAL  n  IN  THE 

NO -STOCKOUTS -ALLOWED  MODEL 


Chapter  II  developed  the  no-stockouts- allowed  model  and 
gave  in  equation  (2-6)  the  minimum  value  of  K  for  a  choice 
of  n,,...,nm.  In  this  chapter  a  method  for  finding  optimal 
n,  ,...»nm  is  given  and,  for  the  case  Aj_  =  A,  for  all  i,  an 
optimal  search  algorithm  is  stated  and  proven. 


A.  THE  GENERAL  CASE 

Let  n  be  the  m-tuple  (n,,...,nm)  of  positive  integers. 
Let  K(n)  be  the  optimal  K  for  each  choice  of  n  as  given  by 
equation  (2-6).  Let  n*  be  the  optimal  value  of  n. 

Minimizing  K(n)  could  involve  considering  a  countably 
infinite  number  of  alternatives  for  n* .  Another  potential 
problem  in  seeking  n*  is  that  there  may  exist  a  sequence 
of  n's,  call  it  In  ,  such  that  the  ]«’*'  Thus, 

—  t  —  3  j  -»oO 

even  though  K(n)  is  positive  and,  hence,  bounded  below,  it 
may  not  have  a  minimum.'*' 

Fortunately,  neither  of  the  above  problems  arise  with 
K(n) .  Theorem  3-1  and  its  corollary  .  given  below  show  that 
n*  lies  in  a  finite,  one-sided  neighborhood  G  of  the  point 
(1,...,1)  whose  boundary  can  be  computed  easily.  Hence,  n* 
can  be  found  by  considering  a  finite  number  of  candidate 
points  n. 


As  an  example  of  this  peculiarity,  consider  y(n)  = 
for  n  =  1,2,3,...  Clearly,  y (n)  >  0  for  all  n;  yet  °* 

Hence,  there  does  not  exist  a  finite  n*  such  that  y(n*)  £  y(n) 
for  all  n. 
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Theorem  3 


- 1 :  For  each  j  5  ni  ,  there  exists  an  Nj 


such  that  K(n)  >  K(l,...,l)  whenever  n^  ^  regardless  of 
the  choice  of  n^,  i  ^  j  . 

Proof :  Let 

=  51  {  KCo)-  ClCif 

Note  that  L.  (i,  ••• ,  i  )  if  and  only  if 

KCO  >  >  K  Cl,-..  1)  .  Let  be  arbitrary 

but  fixed.  Let 


C<  r  Ao  +  I 

X o c o «)D i  ^ X i  D (l 


and 


h  n- 

Then  &  >  O  and  ^6,  >  O  .  Then 

=  (VD*,U  *  i  L  D- 

i  f  n  Xo^Cn—  1)  I-S1>j  n  a- l 

+{°<A  +  o<  ~V~  **  *?•  1 Jj. 


Let 


Mc^j)=  -1)  A^XoDj  + 

Since  LCQ  )  -  M  (\)  >  0  ,  L  (.D  )  >  M  (  *\s )  . 
b  -  Ul,->i)  ,  and  let 

Nr-  r-t[is^5i  +ij^l 


Let 


where  Int  [•]  is  the  greatest  integer  function, 


18 


Since  M(n^)  is  a  strictly  increasing  function  in 

N-,  M(n.)  >  M(N-)  for  n.  >  N . .  But 
J  J  J  3  3 


Nj  > 


b  —  Aj  I  j  Dj 
Ajte  D; 


and 


MIMO  >  )  AitoDj  +  AjXjDj 

=  b  =  L  Cl, •••,!). 


Hence,  for  all  n^  >,  Nj  and  regardless  of  the  value  of 
the  n^  ,  i  #  j  ,  LCQ}  >  M(^)  > 

Then  for  the  same  conditions  on  n,  K(n)  >  K(l,...,l); 
and  the  proof  is  complete. 

Theorem  2-1  shows  that  the  neighborhood 
(jp  =  [  Q  ’•  1  1  A;  <■  Nt  a*<A  1  <  L  <  m  ^  is  a  hyper  box  in 

the  positive  integer  lattice  in  m-space  with  one  corner  at 
the  point  (1,...,1).  The  following  sketch  shows  what  G 
looks  like  if  n  =  (n,,  na)  . 
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An  immediate  implication  of  theorem  3-1  is  the  following 
corollary  that  n*  lies  in  the  bounded  hyper  box  G. 

Corollary  3-2:  For  K  given  by  equation  (2-6) 

Proof :  Either  n*  =  (1,...,1),  or  it  is  not.  If 

n*  =  (1,...,1),  then  n*  =  G.  Assume  now  that 
n*  (1,...,1).  Then  K(l,...,l)  >  K(n*) .  Now 
for  all  n  £  G, 

K(n)  >  K(l,.  .  .  ,1)  >  K(n*) 

which  implies  that  n*  is  not  an  element  of  the  complement 
.  of  G.  Hence,  n*  £  G  .  In  either  case  n*  £  G.  This 
completes  the  proof. 

The  previous  theorem  and  corollary  show  that  in  the 
general  case  n*  and,  hence,  the  optimal  policy  (  Qt  } ) 
can  be  found  in  a  finite  number  of  steps  by  an  exhaustive 
search.  The  exact  number  of  steps  required  is  equal  to  the 
number  of  positive  integer  lattice  points  in  G;  namely, 

TT  For  example,  suppose  that  m  =  5  and 

L=  i 

=  11  for  i  =  1,...,5.  Then  the  number  of  iterations 
required  to  find  n*  is  100,000! 

B.  THE  CASE  Ai  =  At  FOR  ALL  i 

The  balance  of  this  chapter  is  devoted  to  developing  an 
optimal  search  algorithm  to  find  n*  for  the  special  case  of 
having  the  same  ordering  cost  at  each  of  the  lower  echelon 
activities;  that  is,  =  At  for  all  i.  The  method  of  search 
•  is  demonstrated  by  the  following  example  in  which  n  =  (n  , ,  nt)  . 
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Suppose  N,  =5  and  =  3.  The  search  begins  by  com¬ 
puting  K(i,l)  [  actually  L(l,l)j  and  saving  (1,1)  as  the 
first  candidate  for  the  optimal  solution.  The  two  adjacent 
lattice  points  in  G,  (1,2)  and  (2,1)  are  compared  to  find 
out  which  one  gives  the  smaller  value  for  K.  Say  that 
K(l,2)  >K(2,1).  Now  K(2,l)  is  compared  to  K(l,l)  to  find 
out  whether  (2,1)  will  become  the  updated  candidate  for  n*. 

If  K(2,l)  <  K(l,l),  then  (2,1)  becomes  the  updated  candidate, 
Whether  the  candidate  changes  or  not,  the  search  continues 
from  the  point  (2,1).  Now  the  two  points  adjacent  to  (2,1) 
for  which  either  n,  or  na  increases  (that  is,  (3,1)  and 
(2,2))  are  compared  in  the  same  manner  as  were  (1,2)  and 
(2,1).  Say  that  K(2,2)<  K(3,l).  Then  K(2,2)  is  compared 
to  the  candidate  for  n*  as  before,  and  either  an  update 
occurs  or  it  does  not.  In  either  case  the  search  continues 
from  the  point  (2,2).  The  search  terminates  when  either 
n,  =  N,  or  =  N^,  and,  at  that  point,  the  coordinates  of 
the  last  candidate  give  n*.  The  diagram  below  shows  how 
the  above  search  began  and  how  it  may  have  continued. 


„+  o*  Path 

Caw  d,  dale  g 

o>o 

Ci,  1) 

U,0 

0,0 

(2,2.) 

U,2) 

(3,2) 

(2,2) 

(3,3)  S+op. 

hM2,2>. 
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The  algorithm,  itself,  will  be  stated  at  the  end  of 
the  chapter.  Some  theoretical  groundwork  is  required  to 


prove  that  this  method  of  search  works. 

In  the  following  development,  let  g (n)  be  a  positive 
increasing  function  in  each  n^  defined  over  the  positive 


integer  lattice  in  Euclidean  m-space.  Let  n^  =  (1, 


,1). 


•  « 

Let  n^_^  =  (n1,...,nm)  be  an  element  in  a  sequence  of  n's 
which  track  the  progress  of  the  search  through  the  positive 
integer  lattice.  Then  let 

and  let  n^  be  such  that 


^3*  = 

Then  the  sequence  j  n  j  |  traces  the  coordinate-wise,  unit 
step-sized  path  of  steepest  ascent  in  g(n). 

Define  h  (Jl)  by 


Thus,  h  is  defined  on  all  positive  integers  and  corresponds 
to  the  value  of  g (n)  along  the  path  of  steepest  ascent. 


Now,  when  =  A j  for  all  i,  K  contains  a  function 
g(n)  for  which  it  is  important  to  find  the  coordinate-wise 
path  of  steepest  ascent.  Without  justifying  its  existence 
at  this  point,  the  particular  g(n)  will  be  defined  and  a 
property  of  its  associated  function  h  will  be  noted. 
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Theorem  5-5:  Let  g(n)  =  2I»  where  gi(0)  =  0 

~  «u-i  4 

and  gi(n^)  =  J^D^XT  J^TT)  for  ni  >  !•  Suppose  that 

>  0  andl?^  >  0  for  all  i.  If  h(  $  )  is  defined  for  g(n  ^  ) 
as  above,  then  h  C J?  D  ^  g  (n)  for  all  n  such  that 
1-- 


Proof : 


Let 


where 


t 'h. 

Let  h  (0)  =  0.  Finding  the  —  element  of  the  sequence 
is  equivalent  to  picking  the  Jl  largest  elements 
from  the  set  and  letting  n^  equal  one  plus  the  number 

of  elements  picked  from  jjX  ^  .  Then  h  ( )  =  g  (n^). 

For  any  other  n  such  that  ,  g  is 

found  by  picking  the  largest  n^  -  1  elements  o^X*L  ,  taking 
their  sum,  and  then  summing  over  i.  Since  the  sum  of  the 
largest  elements  in  X*  must  be  at  least  as  large  as 

A\. -»  ^ 

L  IT  (j—  largest  element  in  )  such  that 

0  “  1 

^  r  (l  »a-,)-rvi,  u  for  any  2.  such  that  X  =  L(n^  -  m; 

and  the  proof  is  complete. 

Consider  I_»  (n)  with  =  A  for  all  i. 

Ll0)=  (/U*  t  -j. 

r>f 

Let  1^=1^-  IQ  and  suppose  for  the  present  that  1^  >  0 

for  all  i.  Let  the  set  J  =  .  Then 

L»C0)  =  (ao+  A,Lnk|[rftZ.Dt  +  XI;Dw-  15^2  j tjvo } . 

The  last  expression  in  the  second  set  of  brackets  is  g  (n) 

in  the  previous  theorem. 


23 


The  following  theorem  completes  the  proof  that 
coordinate-wise  steepest  descent  works  for  L  (n)  and, 
hence ,  for  K  (  n  )  . 


Theorem  3-4: 


Let  n  ^ 


be  defined  as  before  and  let 


W  -  } 

min  LCn)  -  LlO). 


If 


mm 

n 


L  (  n  )  exists ,  then 


Proof :  First  note  that  theorem  3-1  and  its  corollary 

showed  that  IjCO^  exists.  Let  a t  =  A0  ,  a^  =  A,  ,  and 
a3  =  Io  +  .  Then 

-  [a,  +aiZnvJ  j  CLj 

Clearly,  vvW  m  Lcq)  <  .  Suppose 

a  "*  £ 1¥ 

Lin)  <  rvM'n  Lc.n')  .  Since  the  minimum  exists,  there 

r\  y\  iW 

is  an  n  £  W  such  that  .  Let  (in^- m  . 

—  a  IW  *“*  z' 

A 

Then,  from  theorem  3-3,  0  CO  )  £ 

It  follows  that 

L cS i  *  { n .  *  <u  z  -uj  f4* "  s^'i*  (a.* K - 3C5' >| 


This  contradiction  gives  the  result. 

Corollary  3-5:  ™  ‘  *  K(.D>  =  KCM. 

Proof:  See  the  proof  of  theorem  3-1. 

Theorem  3-3  assumed  that  I ^  >  0.  The  following  theorem 
shows  ivhy  this  assumption  is  not  restrictive. 
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1. 


Theorem 


3-6 :  If  Ij  $  0  for  some  j,  then  nj 


*  = 


Proof:  Suppose  that  nj*  >  1.  Let  n*  be  such  that 

LCD*)  *  X,  (.0 )  •  Now 

'f  If  r  XvDl  XjT>i  1 

Lc**)»  a.+  cul^f  *  n. -  J . 


Let  n  (n,*,...,nj  Since  n.  <  n  * , 

J  J 

/v 

(a  ,  +  azl,  n^)  <  (a  ,  +  a^  £  n*  )  .  Since 

If  <  If  a„d  CI0  E  D,  )• 

Then  L  (ii)  <  L  (  n*)  .  This  contradiction  implies  n  j  *  =  1, 
and  the  proof  is  complete. 

For  the  special  case  of  =  Aa  for  all  i  the  preceding 
theorems  have  proven  that  the  following  algorithm  will 
provide  the  optimal  solution. 


Algorithm 

(1)  Let  J  =  l L'  ^-t>0|  .  If  J  =  $ ,  stop;  n^*  =  1  for 

all  i .  Otherwise ,  set  n  =  Compute  L*  =  L  (  n  ) 


and  N.  =  Int  [  L\  -  +  l]  +  l.  Set  n  =  (1,...,1) 

1  Ie  D*c  — 

and  go  to  step  (2) . 


(2) 


Let  j  £  J  be  that  index  such  that 


for  all  i  £  J.  Increase  n-  by  one.  If  n.  =  N . ,  stop; 

J  ;  J  3 

n*  =  n.  Otherwise  go  to  step  (3) . 

A 

(3)  If  L  (n)  <  L* ,  set  L*  =  L  (n)  and  n  =  n.  Otherwise, 
make  no  changes.  In  either  case,  return  to  step  (2). 
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After  n*  is  found,  then  Q^*,  T^*  and  K*  can  be  computed 
as  shown  in  Chapter  II  (see  equations  (2-6)  through  (2-11)). 

Earlier  in  this  chapter  it  was  shown  that  for  an  exhaustive 
search  over  G  for  which  m  =  5  and  =  11  for  i  -  1,...,5,  the 
number  of  iterations  required  would  be  100,000.  If  this 
algorithm  could  be  employed  (that  is,  if  =  A^  for 
i  =  1,...,5)  for  this  case,  then  the  maximum  number  of 
iterations  required  would  be  50,  an  impressive  reduction. 
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IV.  EXAMPLES  OF  THE  NO -STOCKOUTS -ALLOWED  MODEL 


This  chapter  is  dedicated  to  some  examples  which,  it  is 
hoped,  will  enhance  the  understanding  of  the  model  and  will 
aid  the  reader  in  following  the  use  of  the  algorithm 
developed  in  Chapter  III. 

A.  AN  EXAMPLE  WITH  LOWER  ECHELON  SYMMETRY 

A  two-echelon,  three- activity  inventory  system  is 
operating  with  the  following  parameters: 

AQ  =  $500  ;  A1=A1=  $1  ; 

IQ  =  0.1/year;  1^  =  1^  =  0.05/year; 

C  =  $10/unit;  DL  =  Dj_  =  1,000  units/year. 

The  aim  is  to  find  the  optimal  operating  doctrine.  The 
following  sketch  illustrates  the  system. 
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The  algorithm  in  Chapter  III  terminates  at  step  (1) 
since  I  >  I1  =  It .  Equation  (2-10)  gives 

K*(l,l)  =  $21,002  (or  $1,002  in  variable  costs 
alone),  and  equation  (2-5)  results  in 

Q  *  =  1,002  units. 

Equations  (2-11),  (2-12),  (2-13),  and  (2-14)  yield  the 
remainder  of  the  optimal  operating  policy: 

T  *  =  1.002  years; 

T  *  =  T  *  =  T  *  • 

[0  » 

Qi*  =  QiA;  and 

Q  *  =  2,004  units, 
o 

Had  the  inventory  policy  been  derived  using  the  standard 
EOQ  model  for  each  activity,  the  following  would  have  resulted. 
Since  activities  one  and  two  are  identical, 

Q^*  =  Q^*  =.  63.24  units;  and 
Ti*  =  T^*  =  .063  years. 

This  results  in  a  time-average  variable  cost  for  each  of 
the  lower  echelon  activities  of 

Kt*  =  $31. 62/year. 

To  find  the  time-average  cost  for  the  upper  echelon  recall 

that  T^*  =  T2*  and,  hence,  that  the  upper  echelon  must 

supply  [  Q^*  +  QZ*J  every  T^  *  units  of  time.  The  holding 

and  ordering  costs  divided  by  Tq  for  the  upper  echelon  are: 

Aet>.  ¥  ToC 

K"  Z  ~  (4-i) 
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(4-2) 


Since  Qa*  =  Q^*,  equation  (4-1)  becomes 

K  '  ^75?  +  I.CK-OQf. 

Treating  nt  momentarily  as  a  continuous  variable,  it  is  easy 
to  see  that  equation  (4-2)  is  convex  in  n^.  Minimizing  with 
respect  to  nt  gives  two  candidate  integer  solutions  for  the 
minimization  problem,  nt  =  11  and  n^  =  12.  Of  these, 
nA  =  11  gives  the  lower  cost  value  for  the  problem  parame¬ 
ters.  Then  the  total  variable  time-average  cost  in  this 
second  analysis  of  the  inventory  system  is 

.  r.CCn'-DQ*  +  2  K* 

=  $1,413.40. 

Comparing  the  second  analysis  to  the  first,  the  second 
policy  yields  a  cost  which  is  41%  higher  than  that  given  by 
the  optimal  policy  for  the  model  of  Chapter  II.  This  seems 
only  natural  since  the  model  of  Chapter  II  allows  for 
cooperation  among  the  various  activities  of  the  inventory 
system. 

B.  AN  EXAMPLE  WITH  LOWER  ECHELON  ASYMMETRY 

The  following  example  has  some  differences  in  the 
parameters  of  the  lower  echelon  activities  and  the  resulting 
cost  difference  between  models  is  more  dramatic  than  that 
in  the  previous  example. 

Consider  a  two-echelon,  three- act ivity  inventory  system, 
as  in  the  previous  example,  with  the  following  parameters: 
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Aj_  =  $100; 


AQ  =  $500; 


A,  =  $1; 


IQ  =  0.1/year;  Ij  =  I2  =  0.05/year; 

C  -  $10/unit;  =  100  units/year; 

and  =  1000  units/year. 

Again,  since  I  j  =  4  0;  nt*  =  n2*  =  1.  In  this  case 

the  total  variable  time-average  cost  is 

K*  =  $813.08. 


The  optimal  operating  policy  is: 

Qt*  =  147.83  units; 

Q2*  =  1,478.3  units; 

T  *  =  T  *  =  T  *  =  1.478  years;  and 
1  x  o 

Q  *  =  1626.13  units. 
xo 

As  in  the  first  example,  a  policy  can  be  derived  assuming 
independently  operating  EOQ  activities  and  the  variable  costs 
can  be  compared. 

The  EOQ  operating  policies  for  the  lower  echelon  activ¬ 
ities  are: 

Q  £*  =  200  units; 

Qi*  s  63.24  units; 

Tj_  *  =  2  years ;  and 
T2 *  =  .06324  years. 

This  yields  time-average  costs  for  the  two  activities  of 
K =  $100/year,  and 
K2*  =  $31. 62/year. 


30 


There  are  two  alternative  policies  for  the  upper  echelon  to 
employ.  One  is  to  order  for  both  lower  echelon  activities 
on  the  same  cycle.  This  leads  to  the  constraint 
n1T1 *  =  n2T^*.  The  smallest  values  of  n2  and  n^  which  solve 
this  equation  are  n ^  =  1,581  and  nz  =  50,000,  and  they 

result  in  a  policy  with  ridiculously  high  costs  for  the 
upper  echelon.  The  other  policy  is  to  order  for  the  two 
activities  on  independent  cycles,  and  it  is  this  second 
policy  that  is  derived.  The  variable  cost  for  the  upper 
echelon  as  a  function  of  time  is 

=  A0T,'t[vv]  +  AaI*tLn1T1,,l  +  holding  costs  (t)  .  (4-3) 
Dividing  equation  (4-5)  by  t  gives 


Kw- 


Ao  I*t  [■£.*  J  A°  Tnt  L^-J 


+  ave.  l-xsUirv^  c°sO>(t).  (4-4) 


As  t  grows  large,  Int  [-■%  A]  =  J?.(t)  ,  a  large  integer;  and 

Ini  [gjgj  .  1(0  |t(f)  1 

t  -  > 


Similarly, 


t 


JIM***)  "lTl 

L 


^  • 

Average  holding  costs,  by  an  argument  similar  to  the  one 
above,  stabilize  at 

I.C  (Vila*  +  UC  Ou-lLgi.  # 


Then  li™  |<(0-  K  ,  where 
t-»oo 
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K 


A* _ 

n,T,* 


(4-5) 


+  Ao_  +  I.C  IeCK-l)Q* 


"i.TX 


=  K(n1,n2) 


Equation  (4-5)  is  convex  in  n ^  and  n^  if  nt  and  n^  are 
treated  as  continuous  variables.  Minimizing  K  with  respect 
to  n^  and  n^  gives  n^  =  1.58  and  =  15.8.  Four  alternative 
solutions  for  integer  values  of  (n^n-^)  are  (1,15),  (2,15), 
(1,16),  and  (2,16).  Of  these  four,  the  choice  that  offers 
the  smallest  value  of  K  is  n  j_=  2  and  n  ^  =  16.  For  the 
upper  echelon  alone,  then. 


K(2  ,  16)  =  $1,193.45. 

Adding  to  this  the  time-average  cost  for  the  lower  echelon 
activities  gives  a  system  time-average  variable  operating 
cost  of  $1,325.07. 

The  second  operating  policy  yields  a  system  time-average 
cost  which  is  631  higher  than  that  for  the  optimal  policy 
for  the  model  presented  in  Chapter  II.  With  asymmetry  in 
the  problem,  allowing  the  activities  to  cooperate  pays  even 
greater  dividends. 


C.  AN  EXAMPLE  WHICH  EMPLOYS  THE  ALGORITHM  IN  CHAPTER  III 
This  final  example  shows  more  details  of  how  the 
algorithm  of  Chapter  III  works. 

Consider  a  two-echelon,  three-activity  inventory  system, 
as  before,  with  the  following  parameters: 


32 


Ao  =  $5  ; 

I  =  0.05/year; 
C  =  $10/unit; 


A  i  =  A2  -  $5; 

I  1  =  I ^  =  0 . 1/year ; 

=  100  units/year;  and 
D<^  =  1000  units/year. 


For  n4  = 
N  ^  =  Int 

N  2,  ~  Int 
Since  for  i  = 


n2  =  1,  L (n)  =  $1650. 
1650-50 


[ 


25 


+  1]  +  1  ss  66, 


f1650-  500  ,  -i  1  ,  -I  _  fi 
^  250 -  +  1]  +1-6. 


1,  2  , 


=5/2  and 


UDi  


50/2  ; 


Then 

and 


set  n^  =  2.  Now  n  =  (1,2)  and  L(n)  =  $1700. 

Then  n  remains  (1,1)  and  L*  remains  $1650. 

The  table  below  shows  the  progress  of  the  algorithm. 

In  this  example  nj_*  =  n^  =  1.  The  optimal  policy  could 
then  be  computed  using  the  equations  at  the  end  of  Chapter  II. 
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fit 

rta. 

L  CO ) 

L* 

A 

n 

1 

1 

1  (?  5  o 

1  G  5  o 

(i.n 

1 

A 

ITOO 

/  4 5o 

(i,i) 

1 

3 

1^  1*1 

IG  5  0 

(1,0 

1 

4- 

ails' 

16  50 

(1,1) 

4 

£450 

/6  50 

(1,1) 

£ 

5 

£700 

/6  SO 

(1,1) 

5. 

(p 

Stop,  Sinice  nx  =  jsj^  . 
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V.  THE  BACKORDERS -ALLOWED  MODEL 


A.  MODEL  FORMULATION 

This  chapter  extends  the  model  developed  in  Chapter  H 
by  allowing  the  lower  echelons  to  run  out  of  stock  and  to 
backorder  the  shortage  until  the  next  order  is  received. 
For  each  unit  backordered  a  penalty  cost  is  assumed  to  be 
suffered. 

The  sketch  below  shows  the  net  inventory  during  one 
cycle  for  one  lower  echelon  activity.  For  activities  one 
through  m  a  new  variable,  s^,  the  backordered  quantity  for 
one  cycle,  and  a  new  parameter  ir^,  the  shortage  penalty 
cost,  are  introduced. 


(5-1) 


The  total  costs  for  one  cycle  are  given  by  (5-2) 


+  E  (holding  cost  for  lower  echelon  activity  i) 
+  Z  (shortage  penalty  cost  for  lower  echelon 


activity  i) . 


(5-2) 
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The  holding  cost  is  the  product  of  I^C,  the  unit  holding 
cost  per  unit  time,  and  the  area  under  the  on-hand  inventory 
curve  over  n^  periods.  That  is: 

(XX )  x  (Qi-Si)X  -  ^ '«•  "X  J5T  • 


The  shortage  penalty  cost  is  the  product  of  tt^,  the  unit 
penalty  cost,  and  the  total  number  of  backorders  incurred 
over  n^  periods.  That  is,  n^ir^s^.  Thus,  equation  (5-2) 
becomes : 


J  :  /\0  +  Jj  n;  At  +  C  +  ToC  LX 


Cni-lAn; 


4  t  t.  d  +  z 

^  ir  s. 


(5-3) 


Dividing  equation  (5-3)  by  Tq  gives  the  time-average 
cost  for  operating  the  inventory  system. 


K 


n,  Q, 


C  LDi 


+  IX  L 


3  y)'k 


z 


t)L  A  v  £>, 

<a, 


+ 


n,DiQ> 

2n;D, 


Sv.  * 


StSnXj.  I  .  t 

i^rrs:  s  u 


TTy5tn\D, 

n,Q  , 


(5-4) 


Equation  (5-4)  will  be  denoted  by  K(Q^ , s^ , . . . ,sm ,nt , . . . ,n  ) . 

B.  OPTIMAL  OPERATING  DOCTRINE  FOR  GIVEN  n  ,...,n  . 

m 

A  full  statement  of  the  problem  faced  in  minimizing 
time-average  system  cost  would  be: 

minimize  K  SM...  ,5m,  ■'h,--  >,1«) 

subject  to  Q0=  Qi-Di-X  for  l- ™ ’> 

n^i  a«A  ’irttcje**  fo** 

Qy  O  -for  Lsl»  — >1(V1> 

O  i  S  L  i  Qi  for  L-  > 

and  ykT  =  =  T0  -Cor  Ul,-^  ^  Y-U-,™. 
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Equation  (5-4)  presents  K  with  the  equality  constraints 
already  satisfied.  Full  justification  for  this  simplifica¬ 
tion  is  not  presented  here,  but  the  reader  is  invited  to 
consider  what  information  the  Kuhn-Tucker  conditions  would 
give  in  the  case  of  the  equality  constraints  and  how  this 
would  then  be  used  to  simplify  the  other  first  order  condi¬ 
tions.  Consider  only  the  problem  posed  in  minimizing  K 
given  by  equation  (5-4). 

minimize  K(Qt ,s1 » • • • >sm>nl » • • • >nm) 
subj  ect  to  Q  ^  >/  0  ; 

0  *  5i«  -i£§r  and 

n^  ^  1  and  integer. 

Disregarding  the  variables  n^,...,n  for  the  time  being, 
notice  that  the  constraint  set  is  convex.  If  K  itself  is 
convex,  then  the  entire  problem  is  a  convex  programming 
problem  with  a  global  optimum.  The  issue  of  the  convexity 
of  K  will  be  discussed  later. 

A  A  A 

Let  Qi>s^,...,sm  denote  the  values  of  Q1,s1,...,sm 
which  minimize  K  for  a  given  choice  of  nt , . . . ,n  .  Let 
Qt*  jS^*  ,  •  •  •  ,  sm*  jn^*, .  .  .  ,n^  be  the  optimal  values  for  all 
the  variables.  By  arguments  similar  to  those  given  in  the 
proof  of  theorem  2-1  the  following  can  be  shown  to  be  true: 


Theorem  5-1:  K(Qt>st>...,Sm, K(Qr,sr,...,S«, n, 


Now  Qi»s1,...,s  can  be  found  by  the  calculus 


K (Qi  ,  s^  , .  .  .  , s  )  is  a  global  minimum  for  given  n^  , .  .  .  ,n 


m 


37 


if  the  first  partial  derivatives  of  K  with  respect  to 

A 

Q^,si,...,sm  evaluated  at  Q^,s^  , . . . ,sm,nt, . . . ,nm  are  zero 
and  the  Hessian  of  K  is  positive  definite  (so  K  is  convex) 
for  the  specified  n«,...,n  and  for  all  feasible  values  of 
Ql,si  » •  •  •  »sm»  There  are  several  cases,  however,  in  which 
this  approach  to  solving  the  problem  fails  because  either 
the  partial  derivatives  cannot  be  set  equal  to  zero  simul- 

*  A 

taneously  for  some  0  and  some  s^  such  that 

A  y\  A 

Qi  for  each  i,  or  the  Hessian  of  K  is  not  positive 

definite . 

Taking  the  first  partial  derivatives  of  equation  (5-4) 
with  respect  to  Q  ^  and  s^: 


3j<  _  _  AoD,  x  t  c  r  i  mAv.TX 

5Q,  ‘  n.Q,1-  r  •LoU  ^  aAll>,  ~  qx  — 

y  XwC  |0|T)».  _  J _  y  X'».CS;  ntP(  I  y  Tt;5; 

+  &rrj>i  GUX  u  2 .n.Dt  Q*  *-*  to, 


(5-5) 


;  and 


^  „  IvCs-AJ),  ,  t t;  niD. 

5ii  1  -TvC  -  r  *  •  (5-6) 


Setting  equation  (5-6)  equal  to  zero  and  solving  for  s- 
yields : 


J).  (w.Q.XtC  - 

hC  »rL  D, 


(5-7) 


Substituting  this  result  into  equation  (5-5)  and  setting 
(5-5)  equal  to  zero  gives  a  quadratic  equation  for  Q 
in  terms  of  the  system  parameters. 


3.  Q'2  +  3a  =  O  where 


(5-8) 
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and 


f  Y  ^  Cr 

* 

J  L  JUvD. 

•  -  \ 

1 

m 

wi.  At  D| 

n, 


IoC 


T  (r)V-1)n»'py 
L  Xn;D, 


4- 


L 

l 


Ulth  Kl~  D, 

an,i^c 

^,X>iTvC 

anil, 
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Depending  on  the  values  of  B  ^  and  B2  the  method  for 
finding  > si > • • • > sm  varies  so  that  a  case  by  case  analysis 
will  have  to  be  made.  Note  that  it  is  possible  that  B^  can 
be  negative,  zero,  or  positive,  depending  on  the 
its  numerator, 

AoE,  y  fil  AtT),  y  UiDt 
^ ^  fi ,  a^tTi  c  • 

Note  also  that  if  n.  =  1  for  all  i  then  B ±  =  0 ; 

B^  <  0.  Six  cases  will  be  considered  below: 


CASE 

Bt 

Bz 

ONE 

=  0 

=  0 

TWO 

<  0 

>  0 

THREE 

=  0 

>  0 

FOUR 

<  0 

=  0 

FIVE 

<  0 

<  0 

SIX 

=  0 

<  0 

sign  of 


(5-9) 

otherwise , 


In  the  analysis  of  the  cases,  the  following  property 
of  K  will  be  valuable.  For  a  fixed,  positive  Q*.,  K  is 
positive  definite  everywhere  for  s^_i...,s  ..  If  some 

of  the  s^  =  0  and  Q^>  0,  then  K  is  positive  definite  every¬ 

where  in  the  remaining  s-  >  0.  To  see  that  this  is  true 
„  °  i 

consider  5^55:. 
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 IiC 


for  i  =  j ;  and 
for  i  ^  j  . 


=  o 

Hence,  the  Hessian  of  K  in  any  non-void  subset  of  the  s^'s 
is  a  diagonal  square  matrix  with  all  of  the  diagonal  terms 
positive.  This  implies  its  positive  definiteness  and  "the 

resulting  convexity  in  the  s^'s.  Additionally,  when  the 

3  K 

value  of  s^  that  renders  ^3:  =  0  for  each  i  is  substituted 

into  K  given  by  equation  (5-4),  K  is  then  convex  in  Q^. 

This  method  for  finding  the  minimizing  values  of 
,sA, .  .  .  ,  and  s  is  justified  in  Section  D. 

CASE  ONE :  B ^  =  0  and  B ^  =  0. 

In  this  case  equation  (5-8)  cannot  be  used  to  solve 
for  since  it  yields  no  information.  Let  >  0  be 

arbitrarily  large  but  fixed.  Then  s^  is  given  by  equation 
(5-7),  and  the  positive  definiteness  of  K  in  s^,...,s 

A  A 

implies  that  s^,...,s  provide  a  global  minimum  for  the 
given  Q^.  Now  consider  K  as  a  function  of  at  the  value 

A  A 

s4  , . .  .  ,s  .  Its  value  is 

1  m 

K  s  C  L  Du  Ii  • 

A 

Since  K  is  independent  of  Q^,  could  be  any  positive  value. 
Once  it  is  selected,  s^  could  be  computed  using  equation 
(5-7). 

A 

A  problem  arises,  however,  if  is  chosen  such  that 
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oil  ^  Ti  C 

since  equation  (5-7)  would  give  <  0  for  some  i.  The 
problem  is  resolved  by  considering  the  following  two  state¬ 
ments  of  the  problem. 

A 

Minimize  K  CQi.s . ,s  )  , 

r\  j  vv\O.X  TU;Et 

given  Q,  <  lsV>Sw  y~C  . 

Minimize  K(Qt,st , . . . ,sm) , 

given  Qt  <  ltlJm  fTc"  J  and 


s^  •%,  0  for  all  i. 

Since  the  second  statement  is  a  restriction  of  the  first, 
the  best  answer  to  the  second  can  be  no  better  than  the  best 


answer  to  the  first.  This  implies  that  K(Q^,s^ 


,s  )  for 
’  nr 


*  Wy  o.  X  TCy.fc, 

VI  *  liisw  T-,C 


is  greater  than  or  equal  to  ClDi+ 


and  that  there  is  no  advantage  to  choosing  Q|  that  small. 

^  A  ^y|  ^  y  I  L  y  ^  | 

Thus,  Q,  can  be  picked  as  long  as  Qt  %  i<i<m  j.  c  > 

K 

and  s^  is  then  computed  using  equation  (5-7). 

CASE  TWO:  B ^  <  0  and  >  0. 

A 

In  this  case  can  be  computed  using  equation  (5-8); 


namely  Qi 


■  (•  tel’ 


*via.x 


Provided  that  i‘Lim  n,I-kC  ") 


s^  will  be  non-negative  for  all  i;  and  the  problem  will  be 

^  Max 

solved.  If  Qt  <  1<l4J*  rtX.c 


,  then  si  <  0  for  some  i 


would  result.  Thus,  for  Q<  <  „vV?<vx  z- 

L  It  i  £«»v  n,  It  d 

of  solution  is  invalid. 


,  this  method 
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w>a  x 


trL  ILy  l>, 

Suppose  that  h,rfcC- 

A  A 

for  s.  ,  ...,s  that  minimize  K  results  in 
i*  7  7  m 


and  fixed.  Solving 


-  |  expression  +  C  I  Dl  +  ZjTC.lI>;.+  Qt|reC  Z  ^  ~ 


}<(  i  ~=Z~  <  expression  i 

Q,  l 

a  m  a,  x  K>  - 

Additionally,  suppose  that  Qj  <  i<Um  • 


Then  l<  (Q„  <  minimum  K  subject  to 

/V  A 

Qi  =  Qi  .  In  fact,  as  decreases  from  towards  , 
K  decreases.  Let 


rvia  x 


^1  ~  i<i<m  h,X'wC  ^ 

/v 

and  let  Qj_  decrease  from  Qt  to  Q  ^  =  e^.  Minimum  K  is  taken 
on  where  s^  =  0  for  at  least  one  value  of  i.  Let 


Y« 


■[i  :  si  '  °j  > 

=  [  i  :  S;L  >  0 


and 


Now  as  is  decreased,  the  minimization  of  K  given  QA  <  e ^ 
becomes 


minimize  K 
Sl--  L*  y2 


subject  to  Qi  fixed;  and 

s^  =  0  for  i  Z  . 

The  functional  form  of  minimum  K  is: 

\s  -  —•  f  ^e^-‘  +  v  y  niZKjU  jb\.  ? 

K  <*•  l  A  L  2.n,IwC  ( 


C  IS; 


2L  TLV.1>1  + 


L 


(»v-l)n,l>i 

2n;P( 


T  X  cC  n.Dy) 

ft,2**-  ]• 


(5-10) 


The  first  bracketed  expression  remains  positive  and  increases 
as  the  number  of  elements  in  Y^  decreases.  The  second 
bracketed  expression  also  remains  positive  and  increases  as 
the  number  of  elements  in  Y,  increases. 
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Equation  (5-10)  remains  a  valid  expression  for  minimum 
K  in  sj_,...,s  given  fixed  until  Qt  is  decreased  to  the 
point  where  equation  (5-7)  implies  that  another  s^  =  0  for 
some  i.  Then  the  functional  form  of  minimum  K  changes, 
because  the  sets  Y,  and  Y^,  change.  Once  Y^  and  Y^  have 
been  modified  to  restrict  any  of  the  s^'s  from  turning  nega¬ 
tive,  equation  (5-10)  is  again  valid.  This  method  of 
recursively  formulating  K  remains  valid  until  has  been 
decreased  to  that  value  which  requires  Y2 ,  =  X  and 
Y,  =  £l,...,mj.  Then  the  problem  reduces  to  the  one  given 
in  Chapter  II,  s^  =  0  for  all  i. 

As  i  ranges  from  one  to  m,  let  e^ >  e2>  .  .  .  > ep  be  an 
„  rtiKtD, 

ordering  of  the  ratios  ~  •  r-  .  That  is,  let 

j  l  , 

ef  =  maX  (  1  ) 


el  “ 


e 


max 

max 


UiIvC 

(  niTLtD 
[■MvC 


<eF.i  -  i.«  iy 


a  rui 

}• 


A 

If  Qj_  <  e^,  the  optimal  values  of  ,s^  , .  .  .  ,  and  sm  are  those 
which  correspond  to  the  least  value  of  K  from  among  the 
sequence  of  problems  stated  in  the  above  two  paragraphs. 


The  first  problem  in  the  sequence  is  to 

m  i  r»  i  » x  e  |<^ 

Si*  «-<■  Ya 

Subject  to  O  -Vox-  c£.  Y,  ?  a wd 

i  Qt  ^ 
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where  Yj.  and  Y^  are  the  first  YL  and  Y^  formed. 
For  this  problem 

K  =  |  ^  Et+  EiQs  ^re 

Ht- 1 *$  +  i—-h-ln'uc  ] ’ 
El=  >  a"d 


I  ecr-t 


L*-Yi 

L-On,T>j  .  r  Xil^l 

IiY,  1  * 


Now  K  is  convex  in  Q*. ,  and,  hence,  minimum  K  corresponds  to 


Qi  = 
Qi  = 


otherwise 


ej  such  that  KCC-,)  -  £•*  K(eO 


This  procedure  is  followed  p-1  times  for  the  p-1  inter¬ 
vals  for  Qi,  and  for  each  iteration  a  value  for  minimum  K  in 
Q  i  is  generated.  The  value  for  K  where  s^  =  0  for  all  i  need 
not  be  considered  since  the  last  of  the  sequence  of  p-1 
problems  examines  K(ep) .  From  among  these  p-1  alternatives, 
pick  the  least  value  for  minimum  K.  This  will  give  a  feasible 

A  A 

value  for  Qj.  from  which  feasible  s^  can  be  computed  for  all 
i  £  of  that  iteration.  Of  course,  s^  =  0  for  all  i £ 
of  that  iteration. 

CASE  THREE:  Bt  =  0  and  >  0 . 

In  this  case  equation  (5-8)  cannot  be  used  to  solve  for 
Q  ^  since  ^  ^  is  not  defined.  Let  Qj_  >  0  be  arbitrarily 
large  but  fixed.  Then  s^  is  given  by  equation  (5-7),  and  the 
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positive  definiteness  of  K  in  st,...,s  implies  that,  for 
fixed  ,  st  , . .  .  ,s  provide  a  global  minimum. 

Now  consider  K  as  a  function  of  Qt  at  the  point 

CS1  >  •  •  •  ,  sm)  . 


K 


=  Jj.  [fxpr€SS.'6rs(5-^)J  +  C  l  Dl  +  . 


(5-11) 


For  Qj_  ^  e^,  equation  (5-11)  is  valid  and  minimum 
K  =  I »*v\  K  -  C  2H>1  +•  Z-TL-JD;.  . 

G?,— >  co 

But  suppose  that  Q^<  e  .  Then  equation  (5-11)  is  no 
longer  valid,  and  an  analysis  similar  to  that  performed  in 
CASE  TWO  is  required.  The  only  modification  is  that  E  j 
reduces  to  E5  =  S  •  The  minimum  K  is  picked  from 


Y, 


A  D 


among  the  p-1  alternatives.  If  that  value  of  K  is  less 
■fhan  C-LD^*- then  Qt  ,st  , .  .  .  ,sm  are  the  corresponding 

A 

values  for  which  that  K  applies.  Otherwise,  Qi  =  co  and 

A 

s^  =  (TO  for  all  i.  This  suggests  it  is  best  not  to  go  into 
business . 


CASE  FOUR;  <  0  and  B^  =  0. 

In  this  case  equation  (5-8)  cannot  be  used  to  compute 

A  A 

since  it  implies  that  =  0  which  in  turn  would  imply 
s^  <  0  for  all  i.  As  before,  minimum  K  for  a  given  large  Q* 

A 

can  be  found  when  s^  is  determined  by  equation  (5-7). 

Then  K  =  CLDl  +  ZTLtt);.  +  Qt (ioC L \  •  Since  K  is 

strictly  increasing  in  Q ,  K  is  minimized  as  Q  i  is  minimized. 
But  this  argument  is  only  valid  for  >  e^ .  Hence  one 
candidate  for  minimum  K  is 


(5-12) 


_  _  C _ -  (ni-On.Dl 

K  =  C  iDl  +  S,  [roC  L  2  rr„T>, 


For  Qj_  <  =  0  for  some  i;  and  the  formula  for  minimum 

K  changes. 

Since  s^  =  0  for  some  i,  B  increases;  and,  hence, 

B^>  0.  Now  CASE  TWO  applies.  The  minimum  K  under  CASE 

TWO  conditions  is  computed  and  compared  to  equation  (5-12). 

If  K  in  (5-12)  is  less  than  the  CASE  TWO  solution,  =  e t 
* 

and  s^  is  computed  with  equation  (5-7).  Otherwise,  the 
CASE  TWO  solution  applies. 


CASE  FIVE 


B^  <  0  and  B2  <  0. 


In  this  case  equation  (5-8)  cannot  be  used  to  compute 

f  3*  I  '/* 

Q  since  r  is  imaginary.  As  before,  for  a  given 

large  value  of  Qi  the  minimum  value  of  K  is  given  by: 

K-.  ^  ^expression  (s-^j  +  C ZlX  +  ££'«.  Q«|roCZ ^ ^  j 

(5-13) 

Since  expression  (5-9)  is  negative,  K  is  minimized  for  min¬ 
imum  Qi .  But  equation  (5-13)  is  only  valid  for  *  e^ . 
Hence,  a  candidate  for  minimum  K  is  the  K  given  by  equation 
(5-13)  evaluated  at  =  e^_. 

F*or  <  et  construct  Yj_  and  as  was  done  in  CASE  TWO. 
Here,  K  is  given  by 

i  t  A.t>,  .  j  _  r  ]  +  c  j  tv 

K  "  Q,  I  ”  1  k  lY"Z'  J 


_  _  ^  ^ i- 00:  v  ItCn.Dll 

+  +  QJl.ttln-.S,  +  J. 


(5-14) 
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Now,  if 


remains  negative,  then  the  next  candidate  for  minimum  K  is 
K  given  by  expression  (5-14)  evaluated  at  e^.  This  analysis 
continues  (evaluation  at  e^  would  be  next)  as  long  as 
epxression  (5-15)  remains  non-positive.  If  expression  (5-15) 
turns  positive  at  some  iteration,  then  the  procedures  of 
CASE  TWO  apply  for  the  remainder  of  the  analysis.  Then 
minimum  K  is  picked  from  among  the  p  alternatives,  and 

^  \  A 

Qt  »s^ » •  • • »  and  sm  are  the  values  that  pertain  to  the  K  thus 
picked . 

Although  this  search  is  sufficient  to  find  minimum  K 
it  can  be  streamlined  somewhat  so  that  fewer  alternatives 
need  be  considered  for  minimum  K.  As  K  is  transformed  from 
iteration  to  iteration  and  indices  migrate  from  to  , 

K  is  continuous  in  Q^.  That  means  that  K  for  two  successive 
iterations  when  evaluated  at  the  common  boundary  point  is 
the  same.  While  expression  (5-9)  remains  negative  we  have 
seen  that  minimum  K  in  is  taken  on  for  Q^.  equal  to  the 
left  hand  or  smaller  boundary  value.  Thus,  as  K  changes 
from  iteration  to  iteration  and  as  long  as  expression  (5-9) 
remains  negative,  K  continues  to  decrease  with  decreasing 
Qi.  Accordingly,  past  values  of  K  may  be  discarded  as 
potential  minimum  K  since  they  are  greater  than  the  present 
K. 

The  first  value  of  K  that  should  be  retained  is  the  last 
value  of  K  for  which  expression  (5-9)  is  negative  which 
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equals  K  evaluated  at  the  right  hand  or  upper  boundary  value 
for  for  which  iteration  expression  (5-9)  is  positive  for 
the  first  time. 

CASE  SIX:  B t  =  0  and  <  0. 

This  case  is  the  same  as  CASE  FIVE  except  that  the  term 

is  missing.  If  expression  (5-15) 
turns  positive,  then  this  problem  falls  into  CASE  THREE 
except  that  E^  =  21  -  yv — -  .  The  appropriate 

a  »» 

A  ^ 

values  of  ,s^ , . . .  ,  and  sm  are  found  as  they  were  before. 

C.  THE  SEARCH  FOR  OPTIMAL  nt , . . . ,n 

1  m 

Finding  optimal  nt,...,n  involves  many  calculations. 

The  first  problem  that  must  be  addressed  is,  as  in  the 
no- stockouts- allowed  case,  showing  that  only  finitely  many 
points  need  to  be  considered.  In  general  this  cannot  be 
done  since,  without  some  restriction  on  the  parameters,  as 
the  values  of  n^,...,n  range,  the  solutions  for 

A  ^  A 

Q •  ,s.  , .  . .  ,  and s  found  in  Section  B  cannot  be  controlled 
and  the  value  of  minimum  K  cannot  be  predicted.  However, 
if  the  assumption  that  2  ^  for  each  i  is 

made;  the  hyperbox  G  can  be  constructed. 

Consider  the  following  two  problems: 

(1)  minimize  K(q1,s1  , . .  .  ,sm,n1  ,.  . .  ,njn) 
subject  to  n^  Z-  1  integer; 

Q  1  0  ;  and 

O  <  Q,  -for  aM  *,  W 
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(2) 


minimize  K(Qt ,sL  , . . 
subj  ect  to  ^  1 

Q  i,  >/  ®  > 


•  y  a  y  y  •  •  •  y  t  ) 

m  1  m 

integer;  and 

and 

Q,  for  all  i. 


The  first  problem  is  the  minimization  problem  that  CASES  ONE 
through  SIX  address  with  n^,...,n  fixed.  The  second  problem, 
given  that  2/kTlC  £  JT-'l2' has  as  a  solution  s^  and  Qt 
given  by  equations  (5-7)  and  (5-8)  for  a  given  ni,...,nm> 

The  second  problem  is  also  a  relaxation  of  the  first  problem. 
Let  *  (n^ , .  .  .  ,nm)  be  the  optimal  solution  to  the  first 
problem  for  a  given  n^,...,n  .  Let  (nt , . . . ,n  )  be  the 
corresponding  solution  for  the  second  problem.  Then 

X  * 

Kl  -  •••»  « 

The  following  theorem  shows  that  n*  =  (1,...,1). 

Theorem  5-2:  2A*tI{C  for  all  i, 

then  n*  =  (1 , . . . , 1) . 


Proof:  The  following  argument  shows  that,  if 

for  all  i,  then  whenever  n^  2,  regardless  of  the  choice 

of  n^  , 

K*  $  Kt*(l,...l)  <  Kz*(ni  ,...,nm)$ 
where  K*  =  K^* (n^* , . . . »nm*) ;  and,  hence,  n^*  =  1  for  each  i. 

Substituting  Q-l  given  by  equation  (5-8)  and  s^  given  by 
equation  (5-7)  into  equation  (5-4)  for  K  yields  equation 
(5-16),  the  solution  to  the  second  problem  for  given 

ni’  *  ‘  *  ,nm‘ 
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where 


B  3  =  21  +  lj  lC.il>;.  > 


Bs  (n)  = 


B4Cn)  = 


Note  that  if  n.  ^  2  for  some  i  (or  Bj_  <  0,  the  condition 


that  makes  this  solution  for  K  valid),  then  Bj,  B^Cn),  B5(n) , 
B^,  B^ (n) ,  and  Bg (n)  are  all  positive.  Then,  for  K  given 
by  equation  (5-16), 


kAo  )  >  b3  . 


Consider  now  K(l, .  .  .  ,1)  .  Since  B2  >  0  and  Bj_  =  0, 

CASE  THREE  applies.  Recall  that  for  CASE  THREE 

rv>i  n  |<(Q,)  £  1'^  K(Q,)-  3_ 

S,V-,S^  Qr)0°  s*'">s,n  5 

It  follows  that  whenever  n^  :>  2  for  some  i 


1)  *  Kg* (n)  <  Kt*(n) . 


This  completes  the  proof. 

D.  THE  CONVEXITY  OF  K (Qt , s , . . . , s  ) 

Throughout  this  chapter  the  various  efforts  tihat  have 
been  made  to  solve  for  minimum  K  have  disregarded  the  notion 
that  K  may  not  be  convex  in  ,st  , . . .  ,  and  s  .  In  general, 
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it  is  not.  One  way  to  show  convexity  is  to  show  that  the 
Hessian  of  K,  the  matrix  of  second  partial  derivatives,  is 
positive  definite.  The  Hessian  of  K  is 


•  •  • 


aiK 

2>Q? 

3  0,3  s , 

S0.3St  *  *  * 

^Qi^Srvv 

>ZK 

3*K 

y-K 

35  ,^Q, 

3S,1 

as,  »SX 

e 

• 


3s  ‘ 


Now , 


+  C  T  T  Si2  j  TC;S;ru 

Q?  1  L  L  a«,l>i  *  L  n. 


2ZK  _  ^AoD, 

’  n,Q,s 

_  TU  Wyk,  . 

^  .  Tt_C.nO>.  o  .  .  -  , 

'  o.K-.Q.  L2i>  a"d 

-  O  -Pc<“  Wj  , 

Let  (h-.)  be  the  Hessian  of  K.  Since  x ,  *5, 


=  O 


for  i  ^  j,  (h  — )  is  sparse  (if  m  =  9,  then  (h  —  )  has  only 
28,  out  of  100,  non-zero  entries).  However,  the  analysis  of 
(h^j)  for  a  problem  of  arbitrary  size  is  difficult.  The 
case  for  m  =  1  illustrates  this. 
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Suppose  m  =  1.  Then  (h^)  has  only  four  entries  and  is 
positive  definite  provided  that  hJt  >  0  and  ["k,  ki  j  >  0, 
For  Q.>  0,  hM  >  0.  Now  z 

huKxi-h.i  s  0,4 1—  +  rt,1,v-  xd,  3)4  D, 

If  =  2,  this  reduces  to 

A,  (E^jk  *  A.i.e-x,1 

l  n, 

which  is  negative  for  sufficiently  large.  Hence,  (lu.) 
is  not  positive  definite  for  all  feasible  values  of  and 
sl;  and  K  is  not  convex. 

Although  the  convexity  of  K  is  a  sufficient  condition 
to  claim  that  a  feasible  stationary  point  of  K  is  a  global 
minimum,  it  is  not  a  necessary  condition.  The  analysis  in 
Section  B  showed  that  if  K  was  first  minimized  over 
s^ , . . . ,  and  s  (K  is  convex  in  s  ^, . . . ,  and  s  )  and  if  the 


resulting  cost  function  was  then  minimized  over 

(K(Qi  , st  , .  .  .  ,  s  )  is  convex  in  Q* )  ,  a  global  minimum  could 

be  found.  The  following  theorem  formalizes  that  approach. 

Theorem  5-3:  Let 

KCGM*  k«j, ?.,-.*•»)  *  ™;":.7sor 

and  let 

K  fQ  )  =  vvuV»  . 

Q.y 

Then  K(Qt*,s  *,..., sm*)  =  K(Qt)  ;  Qi*  ^  Qi ;  and  s^  =  si(Q1) 
for  each  i. 

Proof :  Since  KCQj,*  ,S£*  , .  .  .  ,s  *)  is  the  global  minimum, 

K(Qt*,s  *,...,s  *)  <  K(Qi). 
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Now, 


K(Qt)  S  K(Q1,s1*,...,sm*), 
and,  hence, 

K(Q4)  K(Q1*,s1*,...,sm*). 

Therefore , 

KCQt)  =  K(Q±*>s 

Finally,  it  follows  that  Q4*  .=  Q*. ;  and  s-t*  =  s^(Q^)  for 
each  i.  This  completes  the  proof. 

This  method  of  finding  the  optimal  solution  is  used 
throughout  CASES  ONE  through  SIX.  It  suggests  a  weaker 
sufficiency  condition  for  (Q^* ,  s^_*  , .  .  .  , sm*)  than  convexity 
of  K.  This  condition  is  proposed  in  the  following  corollary 

A 

Corollary  5-4:  Let  s^  be  the  value  of  s^  which  minimizes 

A 

K(Qi , s1 , . . . , s  )  for  a  given  value  of  Qj.  .  Then  K(QA)  = 
K(Qi,s1(Q1),...,sm(Qt))  is  convex. 

E.  CLOSING  REMARKS 

This  chapter  has  offered  the  development  of  a  backorders 
allowed  extension  of  the  model  presented  in  Chapter  II. 
Section  B  derived  a  method  that  can  be  used  to  find  candi¬ 
date  points  for  a  local  minimum  to  K.  Section  C  showed  that 
n^*  =  1  provided  that  2A^I^C  2  for  each  i.  Section  D 

showed  the  solution  technique  discussed  in  Section  B  indeed 
provided  the  global  minimum.  In  all,  the  general  backorders 
problem  has  been  solved  except  for  the  problem  of  bounding 
the  search  in  n  when  2A^I^C  <  D-  for  some  i. 
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VI.  CONCLUSION 


A.  SUMMARY 

This  paper  has  presented  two  deterministic  multi- echelon 
inventory  models.  The  assumptions  of  the  first  model  led  to 
the  development  of  an  extension  of  the  EOQ  (no- stockouts  - 
allowed]  model.  Through  the  presentation  of  a  handful  of 
theorems,  it  was  shown  how  the  optimal  policy  could  be  found 
either  by  exhaustive  search  in  the  completely  general  case 
or  by  a  rapid  steepest  descent  search  in  the  case  where  the 
ordering  cost  is  the  same  for  each  lower  echelon  activity. 

From  the  author's  point  of  view  the  three  most  inter¬ 
esting  aspects  of  this  paper  are  the  assumptions  for  the 
steepest  descent  algorithm,  the  effect  of  1^  =  1^  -  Iq  on 
optimal  n  in  the  no-stockouts-allowed  model,  and  the  general 
method  for  finding  optimal  K(g)  in  the  backorders  model. 

The  first  point  of  interest  is  the  way  in  which  the 
assumption  A^  =  At  for  each  i  affects  the  search  for  optimal 
n  in  the  no-stockouts-allowed  model.  That  this  assumption 
streamlined  the  search  is  of  some  interest  mathematically, 
but  it  is  of  additional  interest  from  an  economic  point  of 
view.  In  military  applications  of  inventory  theory  it  is 
usually  assumed  that  the  cost  of  ordering  or  requisitioning 
is  constant  for  any  one  class  of  activities.  For  example, 
among  Navy  stock  points  inventory  "pull"  item  orders  have 
roughly  the  same  cost.  In  the  case  of  afloat  activities 
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it  was  long  assumed  that  the  cost  of  submitting  and  pro¬ 
cessing  any  requisition  was  $7.  The  exact  figures  are  not 
important  here,  but  the  principle  lends  additional  justifi¬ 
cation  for  the  assumption. 

The  second  point  of  interest  is  that  the  observation 
that  1^  =  I  .-I0<  0  implies  that  n^*  =  1.  This  says  that  the 
cost  to  the  system  of  holding  goods  at  the  upper  echelon  is 
at  least  as  great  as  the  cost  of  holding  them  at  that  par¬ 
ticular  lower  echelon.  In  this  case,  the  solution  shows 
that  the  optimization  strives  to  minimize  the  holdings  of 
the  upper  echelon  because  of  the  relatively  high  "opportunity 
cost"  presented.  This  notion  seems  to  support  military 
decisions  to  provide  central  procurement  activities  with 
an  unencumbered  "open  to  buy,"  because  the  "opportunity 
costs"  for  alternative  decisions  are  high. 

The  third  point  of  interest  is  strictly  a  mathematical 
one.  Even  though  the  cost  function  in  the  backorders  -  allowed 
model  is  not  convex  in  all  the  continuous  variables  simul¬ 
taneously,  theorem  5-3  and  its  corollary  show  a  global 
optimum  can  be  found  by  decomposing  the  cost  function 
minimization  problem  into  a  sequence  of  two  convex-program¬ 
ming  problems.  Thus,  for  this  problem,  a  weaker  sufficiency 
condition  than  the  standard  one  of  convexity  has  been  found. 

B.  RECOMMENDATIONS  FOR  FURTHER  WORK 

Specifically,  there  are  two  points  in  this  paper  that 
should  be  re-examined  and,  perhaps,  extended.  The  author 
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feels  that  an  optimal  search  plan,  such  as  the  algorithm  in 

Chapter  III,  can  be  developed  for  the  general  no- stockouts  - 

allowed  model.  For  the  backorders  allowed  model  it  was 

shown  that  if  2A^I^C  £  TC*  for  all  i  that  n^*  =  1  for 

all  i.  A  bound  and  search  scheme  should  be  developed  for 

x 

the  case  where  2A.I.C<  TL  •  D.  for  some  i. 

11  11 

A  general  area  suggested  for  further  work  is  to  extend 
the  backorders- allowed  model  for  the  case  in  which  the 
backorder-penalty  cost  may  be  time-weighted.  Additionally, 
a  "lost  sales"  model  could  be  developed. 

C.  CLOSING  REMARKS 

In  appraising  the  value  of  this  new  model  it  should  be 
remembered  that  it  represents  an  attempt  in  basically  a  new 
direction.  Simple  results  usually  forerun  more  complex  ones 
which  are  often  more  directly  useful.  However,  just  as  the 
standard  EOQ  model  serves  as  a  good  first  approximation  to 
more  complex  deterministic  or  stochastic  single  activity 
models;  so  perhaps  this  model  may  be  found  similarly  useful 
for  multi-echelon  models. 
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